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Abstract 

We investigate the asymptotic tail behavior of massive scalar fields in 
Schwarzschild background. It is shown that the oscillatory tail of the scalar 
field has the decay rate of t~^/^ at asymptotically late times, and the oscilla- 
tion with the period 2ii/m for the field mass m is modulated by the long-term 
phase shift. These behaviors are qualitatively similar to those found in nearly 
extreme Reissner-Nordstrom background, which are discussed in terms of a 

resonant backscattering due to the space-time curvature. 
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I. INTRODUCTION 



The late-time evolution of various fields outside a collapsing star has important im- 
plications for several aspects of black-hole physics. For example, the no-hair theorem |l| 
introduced by Wheeler in the early 1970's, states that the external field of a black hole 
relaxes to a Kerr-Newman field characterized solely by the black-hole's mass, charge and 
angular-momentum. Thus, it is of interest to explore the dynamical mechanism responsible 
for the relaxation of perturbations fields outside a black hole and to determine the decay 
rate of the various perturbations. In addition, the dynamical mechanism of generating per- 
turbation fields ingoing to a black hole is of interest in relation to the problem of stability 
of the Cauchy horizon 0. 

It was first demonstrated by Price , regarding scalar, gravitational and electromagnetic 
perturbations of Schwarzschild black hole exterior, that the fields die off at late time with 
an inverse power-law tail. It has been shown that for a spherical-harmonic wave mode of 
multiple number /, a or t^^'^^ decay tail {t being the Schwarzschild time coordinate) 

dominates at late time, depending on the initial conditions. The same tail behavior occurs 
also for Reissner-Nordstrom black holes and fields of different spin. Recently the linear 
perturbation analyses and the nonlinear simulations have been done by Gundlach et al. 0] 
and 0, respectively. (See also |^.) In addition, the recent treatment of tails in the spacetime 
of a spinning black hole has been done by Barack and Ori 0. 

Many previous works are mainly concerned with massless fields. However, the evolution 
of massive scalar fields will become important, for example, if one considers Kaluza-Klein 
theories, in which the Fourier modes can behave like massive fields. Further, the recent 
development of the Kaluza-Klein idea, such as Rendall-Sundrum model in string theory, 
strongly motivates us to understand the evolutional features characteristic to the field mass 
in detail. 

The physical mechanism by which late-time tails of massive scalar fields are generated 
may be qualitatively different from that of massless ones. In fact, it has been pointed out 
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that the oscillatory inverse power-law tails 

if) ~ t-'-i sin(mt), (1) 



dominates as the intermediate late-time, if the field mass ra is small [|T0| (see also [|TT|)- It 
is clear from Eq. (|T]) that massive fields decay slower than massless ones, and waves with 
peculiar frequency uj quite close to m can contribute to massive tail, while low frequency 
waves can contribute to massless tail. Though the oscillatory power-law form (|l]) has been 
numerically verified at intermediate late times, mM <^ mt <^ l/(mM)^, it should be noted 
that the intermediate tails are not the final asymptotic behaviors; Another wave pattern can 
dominate at very late times, when it still remains very difficult to determine numerically the 
exact decay rate II^jO • 



In the previous paper P], hereafter referred to as paper I, we have analytically found 
that the transition from the intermediate behavior to the asymptotic one occurs in nearly 
extreme Reissner-Nordstrom background. The oscillatory inverse power-law behavior of the 
dominant asymptotic tail is approximately given by 

~ r^/Sin(mt), (2) 

independently of angular momentum wave number /, and the decay becomes slower than 
the intermediate ones. The origin of such an asymptotic tail is expected to be a resonance 
backscattering due to curvature-induced potential. This will be supported by the relation- 
ship between the field mass and the time scale when the tail dominates; For mM <^ 1, 
the smaller mM is, the later the t~^/^ tail begins to dominate {mt ^ l/(mM)^, where M 
is the black- hole mass). On the other hand, for large field mass mM ^ 1, the larger mM 
is, the later the tail begins to dominate {mt ^ mM). Therefore, the time scale when 
the tail dominates will become minimum at mM ~ 0(1), which means that the most 
effective backscattering occurs for such massive scalar fields with the Compton wave length 
1/m nearly equal to the horizon radius M. 

In this paper, following paper I which was devoted to the analysis for nearly extreme 
Reissner-Nordstrom background, we investigate the asymptotic tails of a massive scalar field 



in Schwarzschild background. The purpose of this paper is to determine analytically the 
decay rate and to confirm that the asymptotic tail with the decay rate of t^^/^ is not peculiar 
to the nearly extreme black hole. Nevertheless, it is sure that the form of effective potential 
varies according to the ratio of the black-hole charge to the mass. Then, we also study the 
difference of the asymptotic tail behavior in Schwarzschild background from that in nearly 
extreme Reissner-Nordstrom. In Sec. II, to analyze time evolution of massive scalar fields, 
we introduce the black-hole Green's function using the spectral decomposition method |]12 



Differently from the case of nearly extreme Reissner-Nordstrom in paper I, it is difficult 
to analyze the field equations for the arbitrary field mass m. In this paper, therefore, we 
consider the cases of small field mass niM -C 1 in Sec. Ill and large one mM ^ 1 in Sec. 
IV, and we compare the results with those found in nearly extreme Reissner-Nordstrom 
background. Section V is devoted to a summary. 

II. GREEN'S FUNCTION ANALYSIS 

A. Massive scalar fields in Schwarzschild geometry 

We consider time evolution of a massive scalar field in Schwarzschild background with 
the black-hole mass M. The metric is 

ds^ = -[l ]dt^+ [1 dr^ + r^dn\ (3) 



and the scalar field (j) with the mass m satisfies the wave equation 

,2\ 



(□-m^)0 = O. (4) 
Resolving the field into spherical harmonics 

l,Tn 

hereafter we omit the index / of ip'- for simplicity, and we obtain a wave equation for each 
multiple moment: 

i^,tt - ^,r.r. +Vij = 0, (6) 
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where r^, is the tortoise coordinate defined by 



dr 



I _ 2M ' 



and the effective potential V is 



1^=1- 



2M 



/(/ + 1) 2M , 



(7) 



(8) 



B. The black-hole Green's function 



The time evolution of a massive scalar field is given by 



for t > 0, where the (retarded) Green's function G{r,r';t) is defined as 

92 



+ V 



(9) 



(10) 



Q-f-2 Qy,2 

The causality condition gives the initial condition that G{r,r';t) = for t < 0. In order to 



find G{r, r'; t) we use the Fourier transform 



G{r,,r:-t)e'^'dt. 



(11) 



The Fourier transform is analytic in the upper half u plane, and the corresponding inversion 
formula is 



r-oo+ic 

G{n,ri;t) = -7^ Gin,K;uj)e—'duj 



(12) 



where c is some positive constant. The usual procedure is to close the contour of integration 
into the lower half of the complex frequency plane shown in Fig. 0. Then, the late-time tail 
behaviors which are our main concern should be given by the integral along the branch cut 
placed along the interval —m < ui < m. 
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FIGURES 




FIG. 1. Integration contours in the complex frequency plane. The original integration contour 
for the Green's function lies above the real frequency axis. We choose the value of w on the dashed 
line to be tz7 = Irol and that on the dotted lines to be tz7 = e^^'"\w\. 



Now the Fourier component of the Green's function G{r^,r'^] uo) can be expressed in terms 
of two linearly independent solutions for the homogeneous equation 



dr 



^+u^ -V ]iJi = Q 1 = 1,2. 



(13) 



The boundary condition for the basic solution ipi is that it describes purely ingoing waves 
crossing the event horizon, i.e., 

V^i-e"*^"*, (14) 
at r^, — oo, while tp2 is required to damp exponentially at spatial infinity, i.e., 

^2^e-^'% (15) 



at r^, — > oo, where zu = \J'm^ — uP-. Because the complex conjugate ^\ is also a solution for 
Eq. (p!3D, can be written by the linear superposition 



^2 = a^i + 



(16) 



and the Wronskian is estimated to be 



(17) 



Using these two solutions, the Green's function can be written by 



1 'ipi{K,uj)iJ2{r^,uj) 



The contribution G'" from the branch cut to the Green's function is reduced to 



1 a ~ 



Then the main task to evaluate is to derive the coefficients a and f5. 



(18) 



(19) 



C. Effective contributions to tail behaviors 



Here we point out the at very late times mt ^ 1 the rapidly oscillating term e~*'^* leads 
to a mutual cancellation between the positive and the negative parts of the integrand in 
Eg ■ (p!9D ■ Then, the effective contribution to the tail behaviors arises from the narrow range 



e = 2Mw < 1, 



(20) 



which was also claimed in iTO|,|T^ . 



III. SMALL-MASS CASE 



Introducing the non-dimensional variable defined as 



X 



2M' 



Eq.(|T3]) is rewritten by 



d'^i/j 1 

dx"^ x{x — 1) dx 



+ 



4M2cj2x2 AM'^m'^x l{l + 1) 



(21) 



fx - 11 



X — 1 



x{x — 1) x'^{x — 1) 



^ = {). 



(22) 
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To study analytically the mode solutions ip, the asymptotic matching between the inner 
and outer solutions was successfully used in paper 1 for the potential V in nearly extreme 
Reissner-Nordstrom background. However, to apply the same method to Eg. (p2D , it is 
necessary to assume the field mass to be very small or very large. In this section we consider 
the small-mass case niM <^ 1. 



A. Mode solutions 



1. The inner region 1 < x <C l/m?M'^ 

For the small-mass case such as mM <^ 1, truncating the terms of the order m'^M^, we 
can approximate Eq.(p2[) by 



d^ij ^ f 1 1\ dtp 



dx^ V a; — 1 X } dx 



+ 



= 0. (23) 



{x — ly x{x — 1) 

Now the solution ipi satisfying the boundary condition (|14D can be written using the hyper- 
geometric function F as follows, 

^1 = x^++{x - l)-2*A^m^ 1 - AiMm- 1 - ^ 

- i)-.a, J(l-«M,n)r(-2ri ^ (a-,., A-,_; 1 + 2,; i | . (24) 

where and are 



p2 



and 



/i=^/(/ + -) -SMW (25) 



= ^ ± ± 2iMm, (26) 



respectively, and we used the linear transformation formulas (15.3.6) of |T^ in the second 
equality of Eq. (p^ . If estimated in the region a; ^ 1, we obtain 
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r(l - 4iMm)r(2/i) ^1 r(l - 4iMm)r(-2/i) 



X 



(27) 



r(ir+_)2 r(ir__)2 

which is necessary for asymptotic matching in the overlap region 1 <^ x <^ l/(mM)2 with 
the outer solutions valid in the region x ^ 1. 



2. The outer region x ^ 1 



For large x, Eq.(p2D is approximated by 



dx"^ 



+ 



SM^m^ - /(/ + 1) 



i) = Q. 



(28) 



X x^ 

Then we can write the solutions using Whittaker's functions. The mode solution il)2 satis- 
fying the boundary condition (pIS]) should be 



i'2 = W^,^{2ex), 



where 



If estimated in the region x ^ 1/e, we obtain 

r(-2/i) 



r(i + /i-a 



(29) 



(30) 



(31) 



3. Matching 

We can match the solution (p7|) valid at 1 ^ x ^ l/m^M^ with the solution (^) valid 
atl^x^l/ein the overlap region and determine the coefficients a and (3 in Eg. (p!9|) . 
Then the ratio is given as follows, 

_iT{-2^fT{l + AiMm) 



^ (\ I ^ 



r(l-/i-a)r(K_+)2 

1 r(2/i)2r(l +4iMm) 
>(i + /x-a)r(K++)2 



X (complex conjugate) ^. 



(32) 
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B. Intermediate tails 



The effective contribution to the integral in Eq. (|T9|) is claimed to be limited to the 



range leu — m| = 0{l/t) or equivalently w = 0{^Jm/t), since the rapidly oscillating term 
^-lujt ^^\^[Q^ leads to a mutual cancellation between the positive and the negative parts of 
the integrand (see |TD,S|). The intermediate tails becomes dominant at the late times in the 



range 



when the integral ([ToD should be estimated under the condition 



a 



0{mMVrrd) < 1. (34) 



As was discussed in paper I, the small value of a represents that the backscattering due to 
the spacetime curvature is not effective at intermediate late times. Then, using the condition 
(0), we obtain 



-(|^|,e)--(|a;|,e-^-6) 



I- 



(2/)!2(2/ + l)H 

\uj\~^m 



mMe^'+\ (35) 



which is identical with Eq. (51) in paper I. Therefore, it is obvious that the same intermediate 
tails as Eq. (|l]) dominate at intermediate late times (^3|) , which was numerically supported 
by 00. 

C. Asymptotic tails 

As discussed in paper I, the intermediate tails can not be an asymptotic behavior, and 
the long-term evolution from the intermediate behavior to the final one should occur. The 
asymptotic tail becomes dominant at very late times such that 

» ^17^' (36) 
when the effective contribution to the integral ([19D arises from the region 
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which means backscattering effect due to curvature-induced potentiai dominates. Then, the 
ratio of a to /9 is approximately given by 



\Lu\—*m, 



(38) 



where 



V = ^ ^1..^ ^(4m'M2)''e-^'^'^ + ^ / H'^m^M^y'^e''"', (39) 

i (A — ) i (K+.j 



and 

r(-2u)2r(l -4iMm), o,,,,,, ,-^„ r(2u)2r(l -4iMm), o,,o, „ ,-^„ 



As was shown in the paper I, the contribution from the Green's function to the asymptotic 
tail part corresponds to the integral along the dashed line in Fig. ^ which is approximated 

by 

G'^{r^,rl/,t) c::^ —^^jji{r^,m)tjji{rl,m) ! e*(^''''"'^*^e*'^''rfu; + (complex conjugate), (41) 

47rm,2 Jdashcd line 

where the phase ips is defined by 

and it remains in the range < <y9s < 2-7?, even if cr becomes very large, since we have 

|ry|2 _ |^|2 = > 0. (43) 

/i 



Because the terms e and e^*'^°" in Eq.(41) are rapidly oscillating at very late times, the 



saddle-point integration allows us to evaluate accurately the asymptotic behaviors. Let us 
introduce the variable a defined by 

l_^\'\mtf/\ (44) 
mJ 
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Then, the oscillation terms e^'^'^'^" in Eq.(pD can be approximately rewritten into the 
form 

in the limits mt —* oo and uj/m ^ 1, by keeping a to be finite. Here we have 

= + (46) 

Qj 

The saddle point is found to exist at 



a = ao = — ^ . (47) 



V2 

At the saddle point the parameter e defined by Eq.(^) is given by 

eo = 2M^m2 - ujI ^ 2 ' {mMfl'' , (48) 

or equivalently 

which is equal to the value obtained in the nearly extreme case in paper I. Let us give the 
Taylor expansion of / near the point a = ao as follows, 

fs{a) ~ /,(ao) + ^(a - ao)Vl'(ao) 

= 3 — ^ +3(a-ao)^ (50) 



V2 ) 

Then, the integral in Eq. f^lD can be approximately estimated to be 



^^(r,,^; t) ^ ^(27r)t(mM)5(mt)-t sin(mt - ?,{2TimMf'^{mtf'^/2 - (^.((Tq) + 37r/4) 
2\/ 3 

x^i{r^,m)'ii)i{r'^,m), (51) 

at very late times mt ^1/ {mMY . The oscillation has the period 27r/m, and is modulated by 
the two types of long-term phase shifts. First, the term (3/2)(27rmM)^/^(mt)^/^ represents 
a monotonously increasing phase shift. Second, the term v's(cro) represents a periodic phase 
shift with the amplitude \(ps\ smaller than 27r, which will be studied in the following. Eq.(pTD 
is similar to Eq.(69) in paper I, though the corresponding terms of the phase shift are omitted 
in Eq.(69) in paper I since we have focused our attention to the decay rate. 
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D. Phase shift 



Now we study the periodic phase-shift effect caused by the term ips^uo). From the 
equation 

dips 



da 



0, 



(52) 



we find the maximum and minimum extremes, denoted by (ps+{cr+) and ips^{o-_) of at 
cr = (j+ and cr = (T_, which given by 



(53) 



and 



TT 



27Ta± = 71 ^ — ± 9 + const, 



where 



sin 6 



and 



V 



Expanding \'~f/ri\ with respect to mM <^ 1, we obtain 



1-6, 



and 



/64\ 2 

— ] {mMf (for/ = 0), 



2^'+47r(/!)^ 



(54) 



(55) 



(56) 



(57) 



(58) 



(2/ + l)(2/)!4 

in Schwarzschild background. The phase-shift term 5 is characteristic to the black hole. For 
example, for nearly extreme Reissner-Nordstrom background [H, the value of 5 is given by 
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^nearly extreme 



TT 



50 
23 

24'+37r(/!)2 



mMf (for/ = 0), 

mMf+^ (for/>0). 



From Eqs. 



(2/ + l)2(2/)!4 
I) and (|54D the oscillation amplitude of ips is given by 

~ 27r - V26, 



(59) 



(60) 



and the time interval between the two extremes is estimated to be 



2nA(T = 2n(a^ - a,) = n - 29 



^86. 



(61) 



These equations mean that the smaller the field mass is, the more rapidly the change of the 
phase from the maxixum to the minimum occurs (see Fig. 0). The phase shift oscillates 
at the time scale of l/(m^M^), which is much longer than the period of 27r/m. So it is a 
long-term effect which modulates the basic oscillation. 




20000 40000 60000 80000 100000 



mt 



FIG. 2. Time evolution of the phase ips — <fs- in the case of mM = 0.01 and / = 0. 



IV. LARGE MASS CASE 
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A. Mode solutions 



In this section we consider large-mass case such as mM ^1. It is enough to consider the 
narrow region (pO]) in the same manner with the small-mass case. Introducing the function 
^ defined as 



r-2M 



X 



X — 1 



(62) 



the mode equation (p2[) can be approximately given by 



+ 



AM^u^x^ AM^m^x 



(x-V 



X — 1 



(63) 



since the other terms become negligibly small all over the region x > 1 in comparison with 
the two terms in the potential. Then we can write solutions for Eq.(^) using Whittaker's 
functions. The solutions ipi and ip2 satisfying the boundary condition ( p!4D and (|^), respec- 
tively, are 



(64) 



and 



^2 = a/^^W^.,p (2e(x-l)) 



(65) 



where 



p = -zW4M2tu2 



From Eqs. (^) and (|65D , ip2 can be rewritten as 



(66) 



J r(-2p) ^ 



r(2p) 



ra+p-a 



(67) 



and the ratio of a to /? in Eq.(^) is 



a r(-2p)r(i + p-a) 

/3 r{2p)ra - p - a) ■ 



(68) 
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B. Asymptotic tail 



As was discussed in paper I, no intermediate tails appear for the scalar field with large 
mass, namely mM ^ 1. At late times given by Eg. (p7D , the ratio of a to /3 is approximately 
written by 



Ce 



\uj\~*m 



Xe 



(69) 



where 



C = r(2p)(2m2M2)-V^ 



(70) 



and 



(71) 



Substituting Eq.(|5I?D into Eq.(|TU|), we obtain the contribution from Green's function to the 
asymptotic tail; 

-^ipi{n,m)^i{r'^,m) I e^^'^'^-^^'^'-^'^V^'c/u; + (complex conjugate), 

47rmZ Jdashcd line 



G'^(r*,r:;t) 



where the phase ipi is defined by 



and it remains in the range < < 27r, even if a becomes very large, since we have 
|(|2 _ |;^|2 = T{2p)T{-2p) {e^'^P - e-^'^") > 0. 



(73) 



(74) 



Because the terms e~*'^*, e'^™'^ and e~2piie Eq.(^) are rapidly oscillating at very late times, 
the saddle-point integration allows us to evaluate accurately the asymptotic behaviors. Just 
like the small-mass case in the previous section, using the parameter a defined as Eq. (PD, 
the oscillation terms e«(27ro--'^*-2pin€) Eq.(^TD can be rewritten into the form 



(75) 
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in the limits mt oo and cu/m 1, by keeping a to be finite. Here we have 

2V2mMa 



, V27rmM , 4mM , 
fi{a) + a^ + -—-^ In 

The saddle point can be approximated by 



(76) 



for asymptotically late times 

mt > mM. (78) 

Note that Eq. (^) is identical with Eq. (^Tf ). Evaluating Eq. ([7^) using the saddle point 
integration, we obtain the asymptotic tail behavior 

G^{r^,r^;t) ~ ^(27r)t(mM)3(mt)-i 
2v 3 

X sin(mt - 3/2(27rmM)2/^(mt)^/^ - 4mMlneo - (/?i(o^o) + 37r/4) 
x?/'i(r*,m)?/5i(r;,,m). (79) 



The additional relation (|78D is necessary for the asymptotic tail (|79| ) to dominate. The 
oscillation has the period mjlv: and is modulated by the two types of phase shifts. First, 
the terms (3/2)(27rmM)^/^(mt)-'^/^ + 4mM In eo represent a monotonous phase shift. Second, 
the term v^i(o"o) represents a periodic phase shift since it remains within < y?; < 27r. 
The decay rate of the asymptotic tail in Schwarzschild background is which remains 

identical with that in nearly extreme Reissner-Nordstrom background shown in the paper I. 

C. Phase shift 

In the same manner with the small mass case, we study the periodic phase-shift effect 
which is caused by the term e*"^' in Eq. (|73D . In the large mass limit mM ^ 1, the ratio of 
X and C is approximated by 
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in the Schwarzschild case. On the other hand, in the nearly extreme Reissner-Nordstrom 
hmit, we obtain the large mass limit from Eqs.(59) and (60) in paper I as follows, 



X 

^ nearly extreme 

The oscillation amplitude of ipi is 



(81) 



(82) 



for which Eqs. (|80D and (^Tj) give different values, though it remains same that Aip vanishes 
in the large-mass limit. 



V. SUMMARY 

In this paper we have investigated asymptotic tail behaviors of massive scalar fields 
in Schwarzschild background. We have shown that the asymptotic oscillatory tail with 
the decay rate of t~^/^ and with the period 27r/m dominates, which is identical with that 
in nearly extreme Reissner-Nordstrom background. In addition, we have noted that the 
relation between the value of niM and the time-scale when the t~^/^ tail begins to dominate 
also holds; The smaller the value of mM is, the later the t~^/^ tail begins to dominate 
{mt ^1/ (mM)^). The larger the value of mM is, the later the t~^/^ tail begins to dominate 
{mt > mM). 



As far as the intermediate late-time behavior is concerned, our result agrees with ||10 
However they claimed " SI perturbation fields decay at late times slower than any power law" 
10| , which disagree with the expressions (|5T|) and (0). We believe that their simulation 



m 



was not carried out until enough late times for the asymptotic tail to dominate or there 



occured numerical errors in their simulation. Also, our result is consistent with |[Tl|, which 
was done numerically (both for a test field and for the fully nonlinear cases). 

In terms of the phase shift we can understand more clearly the tail behavior to be a 
backscattering effect due to space-time curvature. In the small- mass case mM <^ 1, the 
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intermediate tail with an oscillatory power-law dominates at the time scale 1 <^ mt <^ 
l/(mM)^ before the asymptotic tail discussed here dominates P,[10[ (see also At both 



the intermediate and asymptotic time scales the frequency of a wave which contributes to 
tails is very close to m. However, the oscillation of the asymptotic tail is modulated by 
phase shift effects which do not appear in the intermediate tail. There exist two types of the 
phase shift, which are monotonous and periodic with time. The periodic phase shift crucially 
depends on the field mass. The shift angle becomes nearly equal to 2tt in small mass limit 
mM — i> 0, and the shift oscillates at the time scale of l/(m^M^), which is a long-term 
effect if compared with the basic oscillation with the period of 27r/m. In large mass limit 
mM oo, the shift angle become very small. In scattering theory, it is well-known that the 
phase shift of wave is connected with the presence of the scatterer. In this case, of course, 
the scatterer is the effective potential of black-hole spacetime. Therefore it is clear that tails 
of massive scalar fields are generated from backscattering due to space-time curvature. The 
shift angle in Schwarzschild geometry is numerically different from that in nearly extreme 
Reissner-Nordstrom geometry, but the qualitative behavior depending on mM is similar. 

We can conclude that the asymptotic tail with the qualitatively same behavior dominates 
both in Schwarzschild and in nearly extreme Reissner-Nordstrom background. Therefore it is 
conjectured that the oscillatory tail caused by resonance backscattering at asymptotic 
late times is a general feature in arbitrary Reissner-Nordstrom background. It remains in 
a future work to investigate what kind of instability of Cauchy horizons is caused by the 
massive tails. 

Finally we give a comment in comparison with the effects induced by the field mass m 
on vacuum polarization of quantum massive scalar fields in the thermal state. It has been 
found that the amplitude of vacuum polarization is enhanced around mM ~ 0(1) in the 
case of nearly extreme Reissner-Nordstrom geometry, while it is not seen in Schwarzschild 
case [Q. The field-mass induced effect on vacuum polarization becomes clear only in the 
nearly extreme limit. The reason will be explained as follows; Vacuum polarization of 
massive scalar fields is due to both thermal excitation induced by black-hole temperature 
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and mass induced excitation. In nearly extreme case, of which the black-hole temperature 
is nearly zero, the mass induced effect becomes significant because the thermal excitation is 
suppressed. On the other hand, in Schwarzschild case the mass induced effect is just hidden 
because the thermal effect by the black-hole temperature can dominate. 

We expect that resonance behavior due to the mass of a field interacting with a black 
hole may exist in various processes as a generic feature of black-hole geometry. In order to 
confirm this conjecture, it is necessary to investigate various processes of massive fields in 
more general black hole models. This will be one of the interesting subjects in black-hole 
physics. 
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